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Abstract. 

High molecular weight polymer solutions have a powerful tendency to deposit adsorbed layers when 
exposed to even mildly attractive surfaces. The equilibrium properties of these dense interfacial layers 
have been extensively studied theoretically. A large body of experimental evidence, however, indicates 
that non-equilibrium effects are dominant whenever monomer-surface sticking energies are somewhat larger 
than fcT, a common case. Polymer relaxation kinetics within the layer are then severely retarded, leading 
to non-equilibrium layers whose structure and dynamics depend on adsorption kinetics and layer ageing. 
Here we review experimental and theoretical work exploring these non-equilibrium effects, with emphasis 
on recent developments. The discussion addresses the structure and dynamics in non-equilibrium polymer 
layers adsorbed from dilute polymer solutions and from polymer melts and more concentrated solutions. Two 
distinct classes of behaviour arise, depending on whether physisorption or chemisorption is involved. A given 
adsorbed chain belonging to the layer has a certain fraction of its monomers bound to the surface, /, and the 
remainder belonging to loops making bulk excursions. A natural classification scheme for layers adsorbed 
from solution is the distribution of single chain / values, P(f), which may hold the key to quantifying the 
degree of irreversibility in adsorbed polymer layers. Here we calculate P(f) for equilibrium layers; we find its 
form is very different to the theoretical P(f) for non-equilibrium layers which are predicted to have infinitely 
many statistical classes of chain. Experimental measurements of P(f) are compared to these theoretical 
predictions. 
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1. Introduction 

High molecular weight polymers are extremely surface-active molecules. Even a weak interaction between a 
single monomer and a surface tends to be magnified into a powerful attraction or repulsion when many of these 
monomers are linked together to form a polymeric chain. It is a remarkable fact that surfaces contacting even 
extremely dilute polymer solutions can readily develop very dense polymer layers. Technologies revolving 
around the properties of either synthetic or biological polymer layers are many and varied, including adhesion 
[1,2], coating [3], colloid stabilization [4], fiber reinforced thermoplastics [5,6], flocculation processes [7], DNA 
microaarrays [8] and biocompatibilization [9]. Motivated both by fundamental questions and by technology, 
understanding and predicting the structure and formation kinetics of these soft layers is a major concern 
of polymer science [10]. A crucial aspect of experimental studies is that adsorbed polymer dynamics are 
typically extremely slow for long chains: an individual chain has many surface stickers and interacts with 
numerous other chains impeding its motion. Irreversibility and non-equilibrium effects are therefore very 
common. The subject of this review is experimental and theoretical work on these non-equilibrium effects, 
and though we consider adsorption from dilute solutions, semi-dilute solutions and melts our main emphasis 
is neutral homopolymer adsorption from dilute solutions. This is the simplest and most widely studied class. 

Polymer-surface adsorption systems are naturally classified according to the mode of adsorption. 
Roughly, there are two classes: chemisorption and physisorption (see figure The clearest example of 
irreversibility arises in chemisorption (figure ^b)) where the covalent polymer-surface bonds which develop 
are often essentially irreversible on experimental timescales. Monomer sticking free energies, e, have values 
typical of covalent bonds which are one to two orders of magnitude greater than kT. Chemical adsorption is 
employed in various technologies where polymers are attached by chemical reactions to solid surfaces either 
from a polymer melt as in the reinforcement of polymer-polymer or polymer-solid interfaces [2,6, 11, 12], or 
from solution as in colloid stabilization by chemically grafting polymers onto particle surfaces [13-15]. 

What is less obvious is why non-equilibrium effects are so widely observed in physisorbing systems, 
even for rather weak sticking energies. Available experimental evidence suggests that irreversibility effects 
become important as soon as e becomes somewhat larger than kT. For example the experiments by Schneider 
et al. [16, 17] for polymethylmethacrylate (PMMA) adsorption onto oxidized silica via hydrogen bonding in 
dilute CCI4 solutions (e w AkT) show essentially frozen-in adsorbed chain configurations. Large physisorption 
sticking energies (e > kT) originate in hydrogen bonding or other dipolar forces, dispersion forces or 
attractions between charged groups. Metal and silicon-based surfaces are usually oxidized and many polymer 
species form strong hydrogen bonds with the surface oxygen or silanol groups [18,19]. Biopolymers such as 
proteins and DNA attach tenaciously to many surfaces due to their many charged, polar and hydrophobic 
groups [8,9,20]. Since hydrogen bonds, for instance, typically have energies of several kT [21,22] it is apparent 
that strong physical bonds are very common. This suggests that whether physical or chemical bonding is 
involved, for long chains irreversible effects may in practice be the rule rather than the exception. 




Figure 1. (a) Schematic of physisorption from a fluid polymer phase onto a surface. Adsorbed chains 
consist of loops, tails and sequences of bound monomers ("trains"). When non-equilibrium effects become 
important, layer structure depends on kinetics of adsorption. This review addresses phyisorption from dilute 
solutions in sections 2, 3 and 4 and physisorption from melts in section 6. (b) As in (a) but for chemisorption. 
In this case chains carry reactive groups which can form covalent bonds (shown black) with a functionalized 
surface. Chemisorption from dilute solutions is reviewed in section 5 and from concentrated solutions in 
section 6. 

To understand non-equilibrium layers, one must identify how they differ from equilibrium layers. The 
theory of fully equilibrated layers is far more advanced, at both the mean- field [23] and scaling [24-31] level 
of description. A main result of these theories is expressions for the decay of the monomer density profile as 



Non- Equilibrium in Adsorbed Polymer Layers 



3 




Figure 2. The two broad classes of polymer adsorption, physisorption and chemisorption, have very different 
values of the parameter Q, the local monomer-surface association rate. Q can be thought of as the conditional 
monomer-surface sticking probablity per unit time, given the unattached monomer contacts the surface. 
Though many systems are in practice mixtures of chemisorption and physisorption, a simplified view of 
monomer free energy as a function of distance between monomer and surface is shown, (a) For physisorbing 
polymers, the activation barrier is very small and and monomer-surface association is very likely upon contact, 
i. e. Qt a is of order unity, where t a the monomer relaxation time. When the sticking energy e exceeds a few 
kT, experiment indicates that chains need large amounts of time to escape the surface, presumably due to 
complex many-chain effects, (b) Chemisorption typically involves a large activation barrier, u 3> kT. Many 
monomer-surface collisions are needed to traverse this barrier, Qt a <C 1. The adsorbed state is also usually 
strongly favored, e ^> kT. 



a function of the distance z from the surface. For adsorption from dilute solutions for example, in the scaling 
picture originally developed by de Gennes [24,25], Eisenriegler et al. [26,27], and de Gennes and Pincus [28], 
each adsorbed chain has sequences of surface-bound monomers (trains) interspersed with portions extending 
away from the surface (tails and loops of size s) with distribution Q(s) ~ s~ n / 5 [29-31] leading to a self-similar 
density profile c(z) ~ z~ 4 / 3 . Experimentally, the existence of an extended diffuse layer is well established by 
a large number of neutron scattering [32-37] and neutron reflectivity [38-40] studies. However a universally 
accepted quantitative test of the predicted density profiles has been difficult to achieve, both due to intrinsic 
limitations of the experimental techniques [41] and to the asymptotic nature of many of the theoretical 
results which are valid in the limit of very long chains. Furthermore, for experiments specifically directed 
at equilibrium, ensuring that equilibrium conditions are realised is difficult when the very non-equilibrium 
effects one wishes to avoid are poorly identified. 

Understanding the origin of the observed deviations from equilibrium for weakly adsorbing systems in 
dilute solutions is a major unresolved puzzle in polymer science. At present it is unclear how universal these 
non-equilibrium effects may be. Various effects have been speculated to play important roles. Kremer [42] 
and de Gennes [43] have suggested that if the solution temperature is below the polymer melt glass transition 
temperature, glassy effects may onset at the interface where polymer densities are similar to those of melts. 
Interactions with the surface might further enhance glassiness. Chakraborty and coworkers [44-47] suggested 
flattened-down chains experience strong kinetic barriers due to local steric constraints which drastically slow 
down dynamics. Ideas related to slow cooperative motions, mutual pinning, development of entanglements at 
the interface and crystalization have been proposed by Johner and Semenov [48] , Sommer [49] , Granick [50] 
and Raviv et el. [51] in a series of comments following a recent theoretical work by Raviv et al. [52] which 
interpreted past experiments [53,54] exhibiting non-equilibrium effects. 

In this review we do not attempt an exhaustive review of the vast body of past research work involving 
strongly physisorbing or chemisorbing polymers at interfaces. Instead, with fundamental issues in mind, our 
aim is to (i) assemble and classify a number of theoretical ideas and numerical simulations which articulate 
the community's current level of understanding of equilibrium and non-equilibrium polymer adsorption, and 
(ii) summarize a number of experimental results which we believe are particularly pertinent and which still 
demand a theoretical explanation. The emphasis is on the simplest case: adsorption of neutral flexible 
polymers from dilute bulk solutions. We also review work on irreversible adsorption from semi-dilute bulk 
solutions and melts, motivated by ideas initiated by Guiselin [55]. Polyelectrolyte solutions, polymers with 
complex architectures and non-flat surfaces are beyond the scope of the present review. 

Physisorption and chemisorption will be carefully distinguished. These are characterized by very different 
values of the local monomer-surface association rate, Q (see figureB). In physisorption, monomer attachment 
is usually essentially diffusion-limited, Q — C/t a , where t a is monomer relaxation time and C is a system- 
dependent constant of order unity [56] . Chemisorption is normally much slower [56-58] with Q values typically 
8 or more orders of magnitude smaller than those of physisorption. The origin of this difference is that 
chemical bond formation usually involves a large activation barrier (see figure [21 . Similarly, desorption rates 
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after chemisorption are usually very small and can be ignored. The effect of dcsorption on physisorbing 
systems is more subtle and is discussed in section 4. The above two classes naturally lead to very different 
adsorption kinetics. This is analogous to bulk polymer-polymer reaction kinetics where depending on Q, 
polymer length N, and solvent conditions, the kinetics are described by one of a range of "diffusion-controlled" 
and "mean-field" kinetic regimes [58-61]. Such regimes also arise for end-adsorbing polymers [56,57,62-67]. 

In section 2 we briefly review the equilibrium picture for dilute solutions and in section 3 we discuss 
experimental evidence for non-equilibrium departures from this picture. Theoretical work related to 
physisorbing non-equilibrium layers from dilute solution is reviewed in section 4. We move to chemisorption, 
again from dilute solution, in section 5. Section 6 addresses irreversibility effects involving melts and semi- 
dilute solutions. We conclude with a brief discussion of the experimental and theoretical outlook. 

2. Adsorption from Dilute Solutions: The Equilibrium Picture 

2.1. Structure of Equilibrium Layers 

This section briefly outlines some central results of equilibrium theories of adsorbed polymer layers. Much 
more extensive reviews can be found in refs. [41,68-71]. In the scaling picture developed mainly in the 
1980's [24-28], each adsorbed chain consists of surface-bound monomers and large loops and tails generating 
a monomer density profile c(z) as a function of distance from the surface, z. Eisenriegler et al. [26-28] showed 
that very close to the surface, in the "proximal" region, the density decays as a power law, c(z) ~ z~ m , where 
the critical exponent m ~ 1/3 represents competion between surface-sticking energy gain, chain entropy, and 
excluded volume interactions. The proximal region decay law crosses over to de Gennes' "self-similar grid" 
regime [24,25,28], c(z) ~ z -4 / 3 , beyond a certain distance h ploK . For z > h pmx the polymer layer can be 
thought of as a semi-dilute solution with continously varying local concentration c(z). In this region the 
correlation length measuring the range of excluded-volume interactions, £ = a~ 5 / 4 c~ 3 / 4 , is proportional to 
the distance from the surface, z, since this is the only relevant length scale: £ « z. Here a is defined to be 
the monomer size. Expressing £ in terms of c leads to t 

(o//ip r ox) (a/z) 1/3 , a < z < h pTOX , _ , T , R _ nN 3/5 fi\ 

(a/z) 4 / 3 , Vo* <z<R F prox 1 ' " ° iV U) 

Unless the bulk polymer concentration, c, is extremely small [72] , then the equilibrium layer height is of order 
the Flory bulk coil radius Rf as indicated in equation ffl. In this same range of c the adsorption isotherm 
exhibits a large plateau, i. e. surface coverage T is weakly dependent on c. 

Even in weakly adsorbing polymer systems, e. g. adosrption through weak van der Waals interactions, the 
value of e is usually of order kT. By studying the adsorption/desorption transition in binary solvent mictures, 
van der Beek et al. [18] estimated the sticking energies per monomer of many typical flexible polymers onto 
silica and alumina surfaces from organic solvents to lie in the range 0.5 to 6kT. Hence the width of the 
proximal region is typically of order the monomer size, h pmx w a, and a clear signature of the proximal 
region is hard to probe experimentally. In the following we consider e of order kT or larger. We remark that 
the net monomer free energy of adsorption e includes both the "stickiness" a monomer feels for the surface, 
but also the entropic disadvantage due to constraining local orientational degrees of freedom upon contact 
with the surface. Thus, crudely speaking one can say the stickiness contribution must exceed a crtical value 
e c representing the entropic disadvantage before a monomer can adorb. Computer simulations show e c is of 
order kT and is lattice-dependent [27]. The real situation is more complex, with various contributions from 
electronic and structural factors such as solvent molecule entropy effects, etc [21]. 

The density decay law of equation reflects a power law distribution of loop and tail sizes. Neglecting 
differences between loops and tails and the details associated with the proximal region, then the loop size 
distribution per surface site is [29-31] 

fi(a) » a~V 11/5 . (2) 

Beyond this, Semenov and Joanny [31] showed that the inner region of the layer, z < z* = aN 1 / 2 , is 
dominated by loops while the outer region, z > z* , is dominated by tails; the resulting density profile obeys a 
z~ 4 / 3 law above and below z*, respectively, but with different numerical prefactors. Support for the scaling 
conclusions of equations (JTJ and is provided by Monte-Carlo simulations of Zajac and Chakrabarti [73], 
de Joannis et al. [74,75], and Cifra [76]. These produce a density decay consistent with the z -4 / 3 law for 
long chains. Zajac and Chakrabarti [73] additionally report agreement with equation (pfl. 

f The cross-over distance /i pr ox and the prefactor in the proximal region density law can be determined by demanding (i) a 
smooth cross-over at /i prox and (ii) the osmotic free energy per unit area, J Rf dzkT/t^ 1 , balances the sticking free energy per 
unit area, eac(a). 
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Complementary to the scaling approach outlined above has been the numerical lattice model of Scheutjens 
and Fleer (SF) [77,78]. This is a self-consistent mean field theory which averages excluded volume interactions 
and thus treats self- avoidance in an approximate manner. This approximation however allows numerical 
solutions for the density profile and for the loop and tail distributions and can additionally describe chains of 
finite length. The mean field description becomes more accurate for solvents near the theta temperature (a 
common case) where self- avoidance is a weak perturbation except for the longest loops and tails. The 
existence of the loop- and tail-dominated regions of the layer was in fact first established by the SF 
model [78,79]. The layer height in the SF theory scales as h ~ A 1/2 [77,78] while the density profile 
decays as c(z) ~ z~ 2 (for sufficiently long chains), different to the z -4 / 3 decay predicted by the scaling 
approach, as shown by van der Linden and Leermakers [80]. Analytical mean field solutions for the density 
profile in the limit of very long chains were derived within the ground-state dominance approximation [81] 
by Jones and Richmond [82]. Going beyond the ground state dominance approximation, Semenov et al. [23] 
subsequently generalized this approach to account for finite length effects to leading order and analytically 
describe the different contributions of loops and tails to the density profile. They found that loops dominate 
for z < zJ^ F = aN 1 / 3 while tails dominate for z > z^ F , similarly to the scaling approach of Semenov and 
Joanny [31]. These new methods have revived interest in analytical and numerical mean field approaches to 
polymer adsorption [83-89]. 

Turning now to experiment, the fundamentals of polymer adsorption at the solid/liquid and air/liquid 
interface have been studied in a vast number of experiments. Research prior to 1993 is reviewed in the book 
by Fleer et al. [41]. Given the strong evidence for nonequilibrium effects (see below), in general one should 
be very cautious when comparing experimental findings to equilibrium expectations. Overall, experiment is 
consistent with the general trend predicted by equilibrium theories regarding structure of the polymer layers 
which were studied, although the fine differences between the mean field and scaling pictures are hard to 
distinguish. Very briefly, measurements of the layer's surface bound monomer fraction as a function of total 
adsorbed amount and molecular weight (MW) by techniques such as NMR [90,91], ESR [92,93], or infrared 
spectroscopy [94] give results fairly consistent with the predictions of the SF theory [41]. The thickness, 
h C xp, of polymer layers has been probed as a function of chain length by hydrodynamic methods [41,95-97], 
ellipsometry [98], and the surface force apparatus [99]. Depending on the method, h cxp is proportional 
to a certain moment of the density profile and many existent measurements are compatible with a power 
law, h cxp ~ N a . Certain studies have favored the SF theory predictions [41,96] while others support the 
scaling predictions [97,99]. For flexible polymer species the total surface coverage T as a function of bulk 
concentration is found to be very weakly dependent on the polymer concentration in the bulk except for very 
dilute solutions, in qualitative agreement with both the scaling and the SF theories [41]. For a given bulk 
concentration, meaurements of T as a function of N in good solvents typically show a weak dependence on 
chain length for large iV [41]. This is consistent with the SF and scaling theories which predict T ~ J dz c(z) 
is dominated by the lower, iV-independent limit. 

Small angle neutron scattering (SANS) and neutron reflectivity have been used to probe the density 
profile. These experiments established the existence of a diffuse extended layer but at present there is no 
general agreement as to the exact form of the density decay. A techical difficulty intrinsic to SANS, as 
dicussed in ref. [41], is its limited sensitivity to the more dilute regions of the layer. Neutron reflectivity 
experiments are also often ambiguous, since multiple density profiles can be fitted to the same data. The 
SANS experiments of Auvray and Cotton [35] and Hone et al. [37] are consistent with the z -4 / 3 scaling 
law. However the results of Hone et al. could also be described by an exponential profile (see also [32]). 
SANS experiments by Cosgrove et al. [33,34] do not support the scaling predictions, but are more consistent 
with the SF theory. Other SANS studies by Rennie et al. [36] are inconsistent with both scaling and SF 
predictions, while neutron reflectivity experiments of Lee et al. [38, 39] and Sun et al. [40] have generated 
data consistent with the scaling predictions. 

2.2. Single Chain Statistics and the Equilibrium Distribution of Bound Fractions 

So far this section has discussed many-chain layer properties. Equally important, and characteristic of 
the layer-forming processes, are properties of individual surface-adsorbed chains in the layer. What is the 
spectrum of configurations of individual chains? According to the scaling picture, a typical chain has ND(s) 
loops of length s or greater, where D(s) = / s °° ds'fl(s') ~ s~ 6 / 5 after using equation Semenov and 
Joanny [100] argue that because of screening effects these are essentially independent blobs and their 2D 
spatial extent parallel to the surface is [ND(s)] 1 ^ 2 as 3 ^ 5 = aN 1 / 2 . This occurs for all scales s; in particular, 
a given chain has of order one loop of length TV 5 / 6 , also of size aN 1 / 2 . Hence a typical chain has a lateral 
size of order aN 1 / 2 , the ideal result (to within logarithmic corrections [100]). 
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Figure 3. Equilibrium probabilty distribution P cq of chain bound fraction, /, in good solvents. For very 
long chains the distribution is sharply peaked at a value / of order unity. For realistic values of TV the 
distribution is rather broad. 



A special role is played by another single chain property, directly reflecting the degree to which individual 
chains are bound to the surface. This is the probability distribution P eq (f) that a chain has a fraction / 
of its monomers touching the surface. This property plays a central role in this review, since its features 
closely reflect whether irreversible effects matter or not. In two independent Monte Carlo studies by Wang 
and Rajagopalan [101] and Zajac and Chakrabarti [73] an equilibrium distribution was found with a single 
peak at a value of / of order unity. To our knowledge, P cq (f) has not been calculated analytically, at least 
at the scaling level. In order to compare equilibrium predictions with experimental measurements of bound 
fractions in non-equilibrium layers in later sections, we have calculated P eq (/) by relating an adsorbed chain 
to a ID unidirectional walk combining an initial tail, a sequence of intermediate steps, and a final tail. The 
result, which is derived in the appendix and is shown in figure^! reads: 

JV-l/8 r 

(/max - /)iV 1/6 



peq (/) 



p[x\ 



» f) 

(a: <C 1) 



(3) 



(/max 

This agrees qualitatively with the numerically measured distribution of refs. [73, 101]. Here p[x] is a scaling 
function which cuts off the distribution at a small-scale-dependent parameter of order unity, / max - The value 
of /max is determined by the small loops with length ~ a, whose number is ~ N per chain. The mean value 
of the distribution, / is located at ~ / max — A" 1 / 6 . The fluctuations around the mean are due to the mass in 
long loops of length A 5 / 6 . There is order one such loop per chain and this diminishes / by an amount of order 
TV -1 / 6 . Equation J3J) has a tail at small / which describes a small population of chains, a fraction of order 
TV -1 / 5 of the total, with / values far removed from the mean: these are the chains with tails [31] of length 
up to order N which determine the layer height w aiV 3 / 5 . Notice in figure we indicate that the distribution 
changes behavior at a minimum value, / m i n , which is also of order unity (to within logarithmic corrections in 
A [100]). The value of / m ; n is determined by the largest loop or tail size in the layer, expressions for which 
have been derived by Semenov and Joanny [100]. We note that chains with / < / m j n values do exist but 
their number is predicted to be much smaller [100]. For a typical value N = 1000, one has A" 1 / 6 re 0.32 
and thus the distribution is in practice rather broad. Broadening effects in theta solvents are expected to be 
even stronger. 



3. Dynamics of Adsorbed Polymers and Evidence for Non-Equilibrium Layers 

The most direct evidence for deviations from the equilibrium picture derives from experiments on the kinetics 
of adsorbed polymers. In this section we review analytical and numerical studies of equilibrium layer kinetics 
and then we proceed to review experiments showing deviations from equilibrium. Particularly revealing 
experiments have been those probing bulk-layer exchange kinetics. 
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Figure 4. Schematic of new chain adsorption in an equilibrium polymer layer (shown as a self-similar grid) 
as described in ref. [102]. (a) Entry: a bulk chain reptates into the layer and makes a first attachment to the 
surface, (b) Spreading: the incoming chain establishes an increasing number of contacts with the surface, 
(c) A typical adsorbed chain configuration is adopted, consisting of trains, loops and tails. A similar picture 
was used in ref. [100]. Chain desorption follows the same path in the reverse order. 



3.1. Theories of Dynamics in Equilibrium Layers 

Compared to static properties, much less is established theoretically about equilibrium dynamics. These have 
been studied for good solvents by de Gennes [43,102-105], Semenov and Joanny [100], and Baschnagel et 
al. [106] for bidisperse solutions. The picture emerging from these works is that the layer of bound chains 
has a certain characteristic equilibration time r cq . This can be thought of as the time after which the chains 
following the distribution P cq (/) of figure |3| are completely reshuffled among themselves. The exchange of 
chains between the bulk and the layer was predicted to be a slower process due to the fact that incoming and 
outgoing chains have to pass through unfavored configurations having a small fraction of bound monomers 
(see figure0J. de Gennes assumed reptation dynamics (i. e. entangled layers) and found the longest relaxation 
time of the layer to scale as [102, 104, 105] 

r cq w t s N 3 (entangled layers). (4) 

Here t s is the relaxation time of an adsorbed monomer which, due to wall-polymer interactions, may be much 
larger than the corresponding time in the bulk, t a [107]. Semenov and Joanny [100] assumed unentangled 
layers and Rouse-Zimm dynamics and obtained r oq « t a N 2 (in their work t s w t a was assumed). 

In equilibrium, layer and bulk chains desorb and adsorb at the same average rate, respectively. In de 
Gennes' picture bulk chains adsorb in two stages (see figure HJ. During the first "entry" stage, the bulk 
chain overcomes the exclude-volume barrier presented by the layer's loops and tails and makes its first 
contact with the surface, i.e. / = 1/N. During a second "spreading" stage the chain continues to make an 
increasing number of surface-contacts, / increases up to f mm , and the chain becomes part of the bound layer. 
When entry is rate-limiting he found that the mean lifetime of an adsorbed chain before its desorption is 
t cx w t a N /(j), where </> is the volume fraction of polymer in the bulk. Semenov and Joanny [100] described 
the dynamics using a similar picture, but assuming unentangled layers and Rouse-Zimm dynamics. They 
obtained a slighlty different chain lifetime, r ox « t a N 2A2 /<j) (to within logarithmic corrections). 

Note that the exchange timescale, t cx , has a weak power law dependence on N rather than exponential 
because the incoming/outgoing barrier is small. The scaling t cx ~ 1/0 reflects the linear dependence on 
concentration of the rate of chain arrival at the surface. Note also that even for the highest dilute solution 
concentrations, = 0*, where 0* = iV~ 4 / 5 is the chain overlap threshold concentration [81], one still has 
r cq <c t gx . A prediction of the above works is that chain desorption into pure solvent, <f> — > 0, is extremely 
slow, which is well-established experimentally [108]. 

Now suppose one starts with a layer of labeled chains in equilibrium and replaces the bulk solution 
with identical but unlabeled chains of the same concentration at t = 0. An important prediction of the 
above theories is that the decay of the surface coverage of labeled chains, T, is a simple exponential for all 
times [43,100]: 

r(t) = r(o) e - t/T - (5) 

An implicit feature of equation JHJ is that there is a single observed desorption rate since r cx 3> r oq , i. e. the 
desporption process is slow enough to sample an average over all equilibrium chain states in the layer. Note 
this result in fact assumes a given desorbed labeled chain does not readsorb, i. e. memory is lost instantly. 
Experimentally, this necessitates a mixing process in the bulk to dilute the desorbed chains. In the absence 
of such mixing diffusion returns a desorbed chain to the surface repeatedly, an effect which can lead to 
non-exponential decay [60,109,110]. 

The kinetics of polymer layer build up starting from empty or "starved" surfaces is more complex and 
has been considered in refs. [43, 100, 102, 111]. 
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3.2. Dynamics of Adsorbed Polymer Layers: Monte Carlo Simulations 

This sub-section provides a brief review of numerical Monte Carlo simulations of dynamics in many-chain 
polymer layers in contact with dilute solutions (for dynamics of single chains see refs. [112-116]). The 
simulations reported here did not include hydrodynamic interactions. The main results of the simulations by 
Wang et al. [101,117] are qualitatively in agreement with the theoretical picture of the previous subsection. 
They found that the lateral dynamics of adsorbed chains up to N = 100 are consistent with Rouse dynamics. 
For sufficiently sticky surfaces (0 < e < 1.5fcT with e c = 0.5/cT) the value of r cx was found to be much larger 
than the lateral relaxation time, even though the scaling dependence on N was the same. This should be 
contrasted with the Semenov and Joanny prediction that the two exponents differ by a small value, 0.42. 
Wang et al. observed non-exponential exchange kinetics arising from readsorption of desorbed chains. 

Lai [113, 118] studied the layer dynamics as a function of e for N < 80 and interestingly found that for 
e ~ lkT (with e c « 0.9kT) the lateral chain dynamics started to slow down and to approach apparently 
glassy dynamics at e w AkT . This result was claimed to be valid despite an omission in the implemented 
algorithm [113,119]. This report is important since it indicates that the value of e is crucial in polymer 
adsorption. 

Zajac and Chakrabarti [120] studied the dynamics for N — 100 and N = 200 and e + e c = 1.8fcT near 
and above = 0*. Their algorithm involved unphysical reptation moves in order to speed up the dynamics. 
In equilibrium they found a distribution of bound fractions similar to the one of figure |3| and observed that 
the internal dynamics of reshuffling of chains between different / values is complex. The timescale for the 
exchange of adsorbed chains by bulk chains was found to be slower than internal equilibration processes. 
Simple exponential exchange kinetics were observed as in equation JSJ. 

Takeuchi [121] also used the Monte Carlo method with extra reptation moves. For surfaces with e l.QkT 
he observed exponential exchange kinetics while for e ~ 0.9kT readsorption effects were important leading to 
non-exponential exchange kinetics. 

3.3. Experiment: Departure from Equilibrium Picture 

The first experimental studies of exchange kinetics on flexible polymers were performed by Pcfferkon et 
al. [122-125] using radioactive labeling techniques. One study [122, 124, 125] involved labeled polyacrylamide 
(PAM) in water adsorbed through hydrogen bonding onto aluminol-grafted glass beads. The beads were 
exposed to a dilute solution of labeled PAM for approximately 30min until time-independent coverage was 
achieved. The labeled solution was then replaced by a dilute unlabeled PAM solution of varying concentration 
c and the amount of labeled PAM which remained adsorbed, T, was monitored as a function of time as shown 
in figure [SJi) . An interesting result of this experiment was that the exchange rate per labeled chain, shown 
in figure |3Jii), was time- dependent and reached a constant value after a cross-over period of ~ 300min which 
was approximately the same for every c. This asymptotic rate was found to increase linearly with c, as 
shown in figure Efiii). The observed spectrum of exchange times disagrees with equation {JJJ and this can be 
interpreted in many ways as follows: 

(i) The observed non-exponential exchange kinetics is a signature of non-equilibrium. Pefferkorn et al. 
[122, 124] argued that the interface is populated with a spectrum of different frozen or slowly changing 
configurations and, consequently, different kinetic properties. (They proposed that the layer consists of 
a flat sublayer of tightly bound chains which exchange slowly, plus a less tightly bound population which 
exchange more rapidly). 

(ii) The layer is not in equilibrium when the exchange experiment starts but it equilibrates after r eq «300min 
which is larger than the layer's preparationg time. The asymptotic exchange rate then becomes constant 
and equal to l/r ox . The fact that asymptotically t cx ~ l/c as seen in figure Efiii) and the fact that 
t cx > 300min as can be seen in figure |5fii), are consistent with this interpretation and the theories 
reviewed in subsection 3.1. Assuming reptation dynamics, equation and given ./V w 1400, this 
implies a relaxation time of adsorbed monomers of order t s « 10 _5 s. This is much larger than monomer 
relaxation times in the bulk, t a « 10~ 10 s. 

(iii) The layer is in fact in equilibrium but its exchange kinetics and internal equilibration processes are much 
more complex than assumed by existent theories, at least for this system. For example, if the equilibrium 
P eq (/) is very broad and chains with different / values have very different exchange times, then the intial 
drop in T will be due mainly to the most rapidly desorbing chains if their desorption times are less than 

7"cq- 

Issues related to surface density of aluminol groups, polydispersity, and effect of pH (this experiment was 
performed at pH=4 where PAM is neutral while many of the surface aluminol groups were positively 
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Figure 5. (i) Time-dependence of surface coverage T of radio-labeled PAM on glass beads in water in the 
presence of polymer solutions at concentration c = 2 X 10 — 3 (a), 5xl0 — 3 (b), 10xl0 — 3 (c), 17xl0 — 3 mg/mL 
(d). Unlabeled curve at top is desorption in the presence of pure solvent. To is initial surface coverage, (ii) 
Specific exchange rate of radio-labeled PAM, from surface to solution, as a function of time for different bulk 
concentrations as in (i). The rate is time-dependent which indicates non-exponential exchange kinetics, (iii) 
Asymptotic apecific exchange rate of radio-labeled PAM as a function of bulk concentration. [Reprinted with 
permission of John Wiley & Sons, Inc. from Pcffcrkorn E., Carroy A., and Varoqui R. 1985 J. Polym. Sci. 
Polym. Phys. Ed. 23 1997, Copyright ©1985 John Wiley & Sons, Inc] 

charged [122, 126]) may further complicate the dynamics. 

In refs. [123, 124] Pefferkorn et al. performed an experiment similar to their PAM studies, but for 
polystyrene (PS) in CCI4 (under good solvent conditions) adsorbing through hydrogen bonds onto glass 
beads bearing surface silanol groups. The measured specific exchange rate of labeled chains shown in 
figure was found to be time-dependent for high surface coverages, but independent of time for low P. 
The interpretation [123, 127] was that for small T the layer consists of flattened-down chains with identical 
exchange kinetics. At higher T an outer layer of loosely bound chains builds up and one returns to the 
phenomenology of the PAM case (figure 0) with time-dependent exchange rates. 

Systematic experiments on PS exchange kinetics were also performed by the group of Granick. In a 
series of experiments, Frantz and Granick [128, 129] studied the exchange of adsorbed protonated PS (hPS) 
by bulk deuterated PS (dPS). Measurements were performed in cyclohexane at the theta temperature of dPS, 
30°C, slightly below the theta temperature of hPS, 34.5°C [129]. The attractive surface was flat oxidized 
silicon. Studying desorption in binary solvent mictures (using a method developed by Cohen Stuart et 
al. [18,130,131]) it was found that [19,132] e w 2kT. After incubating the surface in dilute hPS solution 
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Figure 6. Exchange rate of radio-labeled PS adsorbed on glass beads in CCI4 as a function of 
time under different conditions. Squares: T = 25°C, c = 0.9 X 10 _11 mol/mL, initial surface coverage 
T = 3.4 X 10- 13 mol/cm 2 . Stars: T = 35°C, c = 5 X 10" n mol/mL, T = 4.2 X 10- 13 mol/cm 2 . [Reprinted 
with permission from Pefferkorn E., Haouam A. and Varoqui R. 1989 Macromolecules 22 2677. Copyright 
1989 American Chemical Society.] 



for different agings times, the solution was replaced by dilute dPS solution and the adsorbed amounts of 
hPS and dPS was monitored as a function of time by measuring infrared adsorption spectra (see figure[7Ji))- 
Similarly to the experiments of Pefferkorn et al. [123, 124] the exchange rate was time-dependent. One sees 
in figure |7Ji) that at long times an apparently exponential regime is reached as in figure Efii). The important 
result however is that the time dependence at short times was found to be strongly dependent on the layer's 
aging time. The surface coverage of hPS at short times could be fitted by an exponential law, T ~ e _t / T <= x . 
It was found that t° x reached a constant after aging times of order hours, depending on both N (of order a 
few thousand) and the concentration of the hPS solution (see figure Efii)). This suggests the layer structure 
was contnuously evolving into an equilibrium or metastable state over periods of hours. The large aging time 
value of r° x was found to be exponentially dependent on N for a fixed molecular weight of displacer dPS 
chains. 

The results shown in figs. and are at least consistent with initially non-equilibrium layers exhibiting 
a spectrum of exchange rates, which then age towards equilibrium. A possible model of these events can 
be constructed as follows in terms of the chain-surface contact distriution, P(f). According to theories 
[17,133,134] of the non-equilibrium layer, after adsorption onto a very sticky surface P(f) has a much broader 
distribution than the equilibrium P° q (f). This includes a finite fraction of chains with very small contacts 
and fast bulk-surface exchange kinetics; it also includes a fraction of order one which has the maxiumum 
possible number of contacts, i. e. chains lying fully adsorbed and flattened onto the surface. One would expect 
these latter may have very large r ox values. Thus when desporption is measured, the first r cx values seen are 
those chains with small / values; as time progresses larger and larger / values have had time to desorb. This 
is the r ox spectrum. Only after the layer equilibration time r eq does the distribution of / values reach the 
equilbrium P cq (f) of equation ©. This is much less broad, centered on a mean / value corresponding to a 
r ox value much greater than r cq . Thereafter, a single averaged mean exchange time is observed. It should 
be mentioned, however, that in the experiments of figure \7\ this last phase of reaching equilibrium is not 
manifest; it is possible that this system never reaches a true equiibrium. 

In another study by Schneider and Granick [135], the displacement of hPS by dPS adsorbed on 
oxidized silicon was explored in CCI4 under moderately good solvent conditions. The adsorption energy 
was e k 1.3kT [19,132,135], measured with the method of refs. [18,130,131]. Once again, r° x was dependent 
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Figure 7. (i) Kinetics of exchange of the adsorbed mass of hPS adsorbed on oxidized silicon in cyclohexane, 
normalized by the surface excess just before replacing hPS by dPS solution. Different aging times are shown 
as follows. Squares, 30min; circles, 6h; triangles, 12h; diamonds, 24h. The initial decay of T can be fitted by 
an exponential, T ~ e - */ T cx. (ii) Initial exchange timescale, t<? x , plotted against aging times of the initially 
adsorbed hPS of various molecular weights. A, MW = 1,080,000; B, MW = 575,000; C, MW = 355,000. The 
MW of dPS is 550,000. Squares and circles represent adsorption from solution concentrations c = 1.0 and 
0.1 mg/mL, respectively. [Reprinted with permission from Frantz P. and Granick S. 1994 Macromolecules 
27 2553. Copyright 1994 American Chemical Society] 



on aging times, and it took hours for it to reach what might be a plateau. A main difference with respect to 
the theta solvent case just described above is that at longer times of order hours, the exchange rate became 
essentially zero and no significant exchange could be detected over periods up to lOh (data not shown). Thus 
a fraction of order unity of hPS chains remained trapped in the layers, even when the displacer chains dPS 
chains were much longer than the pre-adsorbed hPS (see figure |S{b)). These findings seem to support a 
picure where the evolution of the layer during aging leads to a state containing more and more kinetically 
trapped chains. This state was reached more rapidly for shorter chains and the trapped fraction depended on 
the length of displacer chains. The authors of ref. [135] argue that it is unlikely that trapping occurs due to 
the slight difference between the adsorption energies of hPS and dPS. We note that Pefferkorn et al. found 
initial exchange rates of order hrs -1 [123] for a far more dilute system involving PS in CCI4 adsorbed onto 
silanol groups (figure EJ). 

We note that experiments by Dijt et al. [136] on "young" PS layers formed in flow from decalin solutions 
onto silica (at theta solvent conditions where e rj 2kT [132, 136]) showed no evidence of significant non- 
equilibrium effects. These workers used optical reflectometry to measure the surface coverage T in mixtures 
of short (MW w 10,000) and long (MW > 100,000) PS chains. Unaged mixed polymer layers with different 
compositions were generated which lead to measurable differences in surface coverage. It was found that 
layers with the same final composition had approximately the same final surface coverage, independently of 
adsorption history. This was not true for transient layers whose coverage however relaxed rather rapidly, in 
times of order minutes. At present it is not clear if and how the apparently fast equilibration of T is related 
to the relaxation of chain configurations probed by exchange experiments using labeled chains which appears 
much slower. 

Another widely studied adsorbing polymer is poly(ethylene oxide) (PEO) for which the glass transition 
temperature (without solvent) is T g ss — 41°C and the melting temperature is T m ss 70°C [137,138]. PEO 
is thus considerably more flexible as compared to PS which has T g ss 90°C [137]. Dijt et al. [136] studied 
mixed layers of long and short PEO chains adsorbed onto silica via hydrogen bonding in aqueous solutions 
using the same method they used for PS and found fast relaxation dynamics for the coverage of unaged 
mixed layers. More recently, Fu and Santore [139, 140] and Mubarekyan and Santore [138, 141] have further 
explored the kinetics of PEO adsorption using TIRF (total internal reflectance fluorescence). Their studies 
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Figure 8. (a) Plot of the initial exchange time constants, Tg X , for the displacement of hPS adsorbed onto 
oxidized silicon in CCL4 by dPS versus aging time. Bottom graph (circles): dPS MW=550,000. Top graph 
(squares): dPS MW=87,000. Open (filled) symbols denote hPS MW=96400 (MW=575,000). The diamond 
indicates a 96400 hPS/550,000 dPS exchange following a 6h aging of a bare surface in pure solvent, (b) 
Trapped fractions of hPS plotted against aging times. Symbols as in (a). [Reprinted with permission from 
Schneider H. M. and Granick S. 1992 Macromolecules 25 5054. Copyright 1992 American Chemical Society.] 



involved adsorption of unlabeled and fluorescently labeled PEO onto silica from water, a good solvent, in 
flow. Microcalorimetry studies [142] indicate e « IkT. It was found [139, 140] that the displacement of 
short PEO chains by long ones was fast for unaged layers, qualitatively consistent with the findings of Dijt 
el al. [136]. However with increasing chain length and aging, the exchange kinetics started to slow down 
and could not be described by a single exponential exchange law as in equation JSJ [138, 140]. For long 
chains a fraction of the layer appeared unexchangeable during the experiment's duration, of order hours (this 
is similar to the trapped fractions observed for PS by Schneider and Granick [135], section |S^b)). It was 
found that after a certain maximum aging time further aging had no observable effect in the kinetics. The 
observed aging effects are summarized in figure El where the initial exchange rate is plotted as a function of 
molecular weight and aging time. Aging in these experiments was performed in a pure solvent rather than in 
contact with the dilute solution. A difference between the findings of the groups of Santore and Granick (for 
PS) is that Mubarekyan and Santore report an initial exchange rate reaching a plateau value more rapidly 
for longer rather than short chains. At long times the kinetics could be fitted by a stretched exponential. 
This may be due to a continuing relaxation of the layer, but also they may be due to the effects of chain 
polydispersity [140]. Polydispersity effects may be very important in these systems due to the tendency of 
long chains to displace the shorter ones, since long chains suffer less entropic penalty when immobilized at 
the surface [106]. In ref. [141] mixed layers of labeled and unlabeled chains were formed and the kinetics of 
fluorescent chain exchange with unfluorescent bulk chains was studied. Provided the layers were aged for a 
period of a few hours, an interesting result was that the kinetics were identical irrespective of the adsorption 
history of the labeled subpopulation. Thus the picture of different classes of chains whose conformation is 
determined by the time of their arrival on the surface appears inapplicable for this system. 

The most direct evidence of irreversible effects has been provided in another series of experiments by 
Granick's group [16,17,143,144] who followed the adsorption of protonated PMMA (hPMMA) and deuterated 
PMMA (dPMMA) onto oxidized silicon via hydrogen bonding from a dilute CCI4 solution (pure PMMA has 
T g w 100°C). The temperature was 30°C, slightly above the theta temperature of hPMMA, 27°C [143]. This 
is a system with somewhat larger sticking energy, e w AkT [143]. These PMMA layers exhibited extremely 
slow exchange rates with bulk chains: for layers of hPMMA aged for rs lh and exposed to dPMMA solution, 
Tg X w lOOh was measured, much longer than the duration of the experiment [143]. Infrared absorption 
spectroscopy was used to measure simultaneously the total coverage T and the surface-bound part Tbound 
as they evolved in time [16,17]. This revealed the functional relationship Pbound(r), which was found to be 
independent of molecular weight, see figure HUT a^. Using deuterium labeling in combination with infrared 
dichroism measurements they could measure the bound fraction / of chain subpopulations. It was found 
that early arriving chains had much higher / values than late arrivers, and these / values were frozen in 
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Figure 9. Values of the initial exchange rate, fc ex ~ 1 /( T ex c )> f° r * ne displacement of labeled PEO chains 
adsorbed onto silica in water by unlabeled chains of the same MW, as a function of aging time and MW. 
The curve through each data set is a single exponential, k ex (r) = fc ex (0)e _ */ Ta s° , summarizing the layer 
evolution for each particular molecular weight. The inset summarizes the aging time constants r ago for the 
single-exponential fits in the main part of the figure. [Reprinted with permission from Mubarekyan E. and 
Santore M. M. 2001 Macromolecules 34 4978. Copyright 2001 American Chemical Society.] 

during the experiment's duration. Measuring the distribution of / values among chains they found a bimodal 
distribution shown in figure EHb),(c) with two peaks at small and large /, respectively. This is strikingly 
different to equilibrium layers where the distribution exhibits a unique peak shown in figure [3J for very large 
N all chains within an equilibrium layer become statistically identical and to within small fluctuations are 
characterized by the same value, / = Fbound/T- The results of other experiments by Soga and Granick [145] 
studying poly(dimcthylsiloxane) (PDMS) adsorbing onto oxidized germanium from CCI4 were consistent with 
the existence of a bimodal distribution. 

Parallel to the above works on exchange kinetics, evidence for very slow and possibly non-equilibrium 
behavior has been gathered in experiments measuring forces between polymer covered surfaces using the 
surface force apparatus. These systems have been modeled using equilibrium concepts [146-149]. Significant 
hysteresis effects are, however, common: for PS adsorption onto mica in near-theta solvents, Hu and 
Granick [150] and Ruths et al. [151] have observed that following the first compression run, the force 
as a function of surface separation during decompression is different and corresponds to layers of smaller 
thicknesses. Ruths et al. observed that force profiles upon recompression did not relax back to the original 
profile for days. These workers assigned the observed phenomenolgy to extremely slow relaxation processes. 
(For these systems the authors estimate [151,152] e ~ 0.3fcT). Similarly slow dynamics were observed by 
Ruths et al. [153] for poly(ethylene-propylene) on lipid coated mica in good solvent. 

For the more flexible PEO systems, hysteresis effects are also observed but appear to relax faster. 
Luckham and Klein [53] and Raviv et al. [54] studied PEO adsorbed onto mica in toluene, a good solvent 
(in ref. [52] it was estimated that e « 0.25fcT after using the theory of ref. [148] to relate osmotic pressure 
to sticking energy). The layers were incubated in solution overnight, the two surfaces brought into contact, 
and the force measured as a function of separation using a surface force balance. The force profile upon 
decompression or rapid subsequent recompression was found to be different to the initial profile, describing 
a layer of a smaller thickness than the one during the first compression. Recovery of the original force 
profile occured more rapidly for shorter polymers. Recovery times for the longer polymers (molecular weight 
310,000) was of order an hour. Theoretical work by Raviv et al. [52] assuming irreversibly frozen chain 
configurations of the compressed layer is consistent with the measured force-distance relationship immediately 
after compression. Johner and Semenov [48] proposed that a dense phase of PEO nucleates on the surface. 
There is also evidence, albeit controversial [154,155], that PEO forms aggregates in solution which adsorb 




Figure 10. (a) Mass bound by hydrogen bonding to the surface, Pboundi versus total mass adsorbed, T, for 
PMMA adsorption onto oxidized silicon from CCU at concentration = O.Olmg/ml. Dashed line: results 
from uninterupted adsorption experiments, MW=107,000. Symbols: results from dosage experiments (see 
refs. [16,17]). Open circles: MW = 90,000, diamonds: MW= 107,000, filled squares: MW=64,000, open 
squares: MW=7700, stars: MW=64,000 at <p = 0.08mg/ml, plus signs: MW=107,000 in toluene solvent, 
filled circle: MW=107,000 from scmidilute solution, (b) Histogram of distribution P(f) of chain bound 
fractions /, derived from uninterrupted adsorption experiments shown in (a), (c) As in (b) but with P(f) 
derived from dosage experiments. [Reprinted with permission from Douglas J. F., Schneider, H. M., Frantz 
P., Lipman R., and Granick S. 1997 J. Phys.: Condens. Matter 9 7699. Copyright 1997 IOP Publishing 
Ltd] 



onto surfaces [156,157]. If so, this could be an important complication in all PEO systems. 

Further information on the dynamics of adsorbed polymers in solution arises from experiments on 
systems which are by their nature far from equilibrium. This includes studies of the kinetic build up of 
layers on initially bare surfaces. A recent review of kinetics in such systems is given by Cohen Stuart and de 
Keizer [108]. Another class of far from equilibrium systems involves the displacement kinetics of adsorbed 
chains by chains of a different species. Experiments performed by the Granick group [158-161] involved 
PS adsorbed onto oxidized silica over a period of 1 hour in a good solvent (CCI4), displaced by a more 
strongly adsorbing polymer: Johnson et al. [158, 159] studied displacement by PMMA while Schneider et 
al. [160,161] studied displacement by polyisoprene. These studies found that coverage of surviving PS chains, 
Tps, obeyed long time stretched exponential laws, Tps ~ e~ (t / T where r' was of order hours and strongly 
temperature dependent. Dijt et al. [162] studied displacement of PS by poly(butyl methacrylate) (PBMA) 
and poly(tetrahydrofuran) (PTHF), and displacement of PBMA by PTHF. Adorption occured on oxidized 
silicon in decalin solutions. Their results suggested that chain flexibility plays important roles in determining 
the magnitude of exchange timescales. Finally, evidence for slow relaxation processes in adsorbed polymer 
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layers also comes from ESR measurements by Pan et al. [163] for adsorption of PS onto porous silica in CCI4. 
4. Theoretical Work on Non-Equilibrium Physisorption 

In the previous section we reviewed experimental evidence for non-equilibrium effects in polymer layers in 
solution. In summary it appears that, with increasing aging, chains adsorbed on strongly attracting surfaces 
undergo internal relaxation processes and become progressively more attached to the surface and harder to 
remove. This internal relaxation process becomes slower for even more sticky surfaces, and for systems such 
as PMMA on oxidised silicon the chains appear to become effectively frozen into configurations which depend 
on their adsorption history. The universality of the existent observations is unclear. 

In this section we review theoretical work addressing these issues. Despite many proposals, the nature 
of the very slow dynamics is not well understood. 

4-.1. Single Chain Physisorption 

Perhaps the theoretically simplest problem to attack theoretically is single-chain adsorption in the limit 
of complete irreversibility of monomer-surface bonds. For physisorption, each attachment is typically very 
rapid after monomer-surface contact, requiring a small number of monomer-surface collisions. We are not 
aware of first principles analytical theories addressing this problem which account for the complex polymer 
hydrodynamics near an interface, subject to an increasing number of surface-polymer bonding constraints. 
(However, somewhat related issues have been addressed theoretically in the context of chain collapse kinetics 
in a poor solvent [164]). Ignoring hydrodynamics, one might expect an upper bound on the adsorption 
timescale r a( j s to be the bulk coil relaxation time Tbuik [81,165], since the constraints formed by attachments 
would presumably accelerate the rate of arrival of new monomers onto the surface. However, the validity of 
even this statement is no longer clear when polymer-surface hydrodynamics are important. 

Numerical studies have provided important insight. Shaffer [166] used the bond fluctuation model to 
study single chain adsorption on a surface with sticking energy lOfcT. Due to the high cost in computer time 
hydrodynamic interactions were not included. Since there was no activation barrier for monomer adsorption, 
his results correspond physically to almost irreversible physisorption. Using N < 80, he found that in the 
presence of excluded volume interactions (physically corresponding to dilute solution) r ac i s ~ -ZV 1,58 . In the 
absence of both excluded volume interactions and hydrodynamics (which would correspond physically to 
adsorption of a few stick chains from a non-adsorbing polymer melt matrix), r a( js ~ -/V 15 . Note that both 
timescales are weaker powers of N than the corresponding bulk relaxation times [81, 166]. Similar scalings to 
those of Shaffer were independently obtained by Ponomarev et al. [167] who used the bond-fluctuation model 
for N < 100. The kinetics of single chains adsorbing onto strongly attractive surfaces have also been studied 
by molecular dynamics simulations by van Eijk et al. [168] and Liu and Chakrabarti [114] 

Konstadinidis et al. [169] used Monte-Carlo simulations to study the structure of a single fully and 
irreversibly collapsed chain on a surface. In the collapsed state, whose nature was dependent on the adsorption 
kinetics, steric and topological constraints at the monomer level prevented every chain unit from actually 
touching the surface. The fraction of a flattened chain's monomers which are physisorbed to the surface was 
measured to be us m 0.85. This should be compared with the value that this simulation would have produced, 
u> — y 1, if the system were allowed to reach equilibrium by introducing, say, a very small desorption rate. 
This is because in the limit e 3> kT the nature of this particular lattice model allows the chain to completely 
flatten down into a two-dimensional "pancake" [43,81]. 

An important open issue is how to go beyond complete irreversibility and study, for very sticky surfaces, 
the dynamics of a single chain after flattening down onto the surface. These dynamics are relevant to layer 
formation on sticky surfaces during the early stages when the layer consists of isolated single chains. In a 
simple picture, one might expect that for e > kT monomer relaxation rates are locally slowed down by a 
Boltzmann factor ~ e( e+u " kT associated with local kink formation laterally or normal to the surface (here 
u is the height of the adsorption potential barrier, see figurel^b)). However, the scaling with respect to N 
would remain a power law. f On the other hand, dynamics may be more fundamentally modified both by 
local steric and global topological constraints. Due to local constraints the degree to which a monomer in 
this flattened configuration can move laterally or desorb is likely to be strongly dependent on the adsorption 
state of nearby monomers. Thus its motion might necessitate the cooperative motion of its neighbors which 
introduces kinetic barriers. Chakraborty and coworkers [44-47] studied numerically a one-dimensional Ising 

f For weakly adsorbed chains, this scaling has been measured [170—174] and numerically simulated [112, 115]. In ref. [170, 171] 
using poly(ethylene glycol) dilutely adsorbed onto self-assembled hydrophobic monolayers, condensed octadecyltriethoxysilanc, 
with e ps 0.5 — lkT, simple power laws were measured for the chain's self-diffusivity: D ~ TV -3 / 2 . 
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model of the effect of such kinetic barriers. For e > kT stretched exponentials were found in the dynamics 
of correlation functions. Relaxation times increased with increasing e and in certain cases diverged as e 
increased towards a value of order kT. These results are similar to those otained in kinetic Ising models of 
glasses [175]. They were intepreted to indicate the onset of glassiness in strongly adsorbing systems. The 
effect of global topological constraints on the dynamics has not been analyzed to our knowledge. 

4.2. Many-Chain Adsorption 

Although at present there is little theoretical work addressing the precise nature of the extremely slow 
dynamics in adsorbed layers, past work has suggested a number of potentially important factors as follows. 

(i) Adsorbed monomer relaxation times may be considerably lengthened due to Boltzman factors e ( e + u )/ kT 
associated with desorption (see fig. [2J. 

(ii) Topological constraints may be important, as suggested in refs. [48,50,51,138]. It is natural to expect 
that if the equilibrium layer is dominated by loops of order the monomer size, then mutual pinning 
down of chains may become very restrictive. Equilibrium theory does indeed predict that when e reaches 
values ~ kT or greater most of the mass of a typical chain belongs to such small loops. This hints at a 
possible transition to more glass-like behavior at this point. 

(iii) Equilibrium theory predicts that for e > kT, the monomer density at the surface (z ss a) becomes of 
order the density in a polymer melt, see equation QJ. Kremer [42] and de Gennes [43] have suggested 
that if the solution temperature is below T g for the melt (as it would be, for example, for PS and 
PMMA at room temperature), glassy effects might onset at the interface where polymer densities are 
similar to those of melts. Interactions with the surface may further enhance glassiness. Consistent with 
this picture, increases of glass transition temperatures have been observed for thin polymer melt films 
deposited onto strongly adsorbing surfaces [176,177]. 

(iv) de Gennes [178, 179] suggested that a dense polymer phase may form near the sticky surface for water 
soluble polymers such as PEO. This would block interfacial chain motion. Johner and Semenov [48] 
proposed an analogous mechanism for PEO in an organic solvent but this was disputed by Raviv et 
al. [51]. Analogous phenomena might occur for other polymers which tend to form crystalline phases. 
Chain stiffness may also enhance local chain alignment and surface ordering effects. 

(v) The dynamics in thin polymer melt films are often much slower than in the bulk, even above T g [177, 180]. 
The theory of Semenov [181] and of Semenov and Johner [182] for confined polymer films predicts 
relaxation times which depend exponentially on molecular weight, similar to the activated reptation 
model of bulk entangled melts [183-185]. Unlike the unconfined case, Semenov found exponential 
dependence for even short entangled polymers in the confined situation. Semenov and Joanny [100] 
proposed that a similar phenomenology might apply to adsorbed polymers in solution. 

Most theoretical studies of non-equilibrium effects have considered the limit of complete irreversibility 
upon attachment. This is directly relevant to more sticky systems such as PMMA on oxidised silicon [16, 17]. 
Barford et al. [186] studied irreversible physisorption, motivated by experiments on polyelectrolyte adsorption. 
Since the limit of high salt concentration was considered, this analysis is also relevant to neutral polymers. 
The build up of the polymer layer profile c(z, T) was followed as a function of increasing surface coverage 
r. An incremental increase in adsorbed amount dT contributes a change in density dc(z,T) whose form 
was determined by numerical solution of self-consistent mean-field equations in which the excluded-volumc 
repulsion due to earlier arriving chains was assumed to be represented by a fixed external field. Physically, 
this models a situation where new chains are irreversibly frozen in after exploring the phase space which was 
available upon their arrival. The incremental profiles were integrated up to the saturation surface coverage, 
r°°, to obtain the final profile c(z, T°°) = dc(z, V). 

The results of Barford et al. [186] for irreversible adsorption compared to equilibrium self-consistent 
mean- field theory are illustrated in figure ^2 The asymptotic surface coverage was found to be less than the 
equilibrium value (see figure ITTT a) ). Relative to the equilibrium density profile c(z), the irreversible profile 
was larger at small and large z but smaller for intermediate z values (see figure ITTT b')'). Another prediction 
is that the final profiles depend on sticking energies, though results for very high sticking energies should 
be interpreted with caution since the equilibrium self-consistent theory predicts r°° — > 00 for infinite chain- 
surface sticking energy (see figure UlT aH. In a subsequent work by Barford and Ball [187] the condition of 
frozen density profiles was relaxed and only the density at the interface was assumed to remain frozen. The 
final conclusions remained very similar. 

Douglas et al. [17] modelled irreversible adsorption by simulating a random sequential adsorption 
process [188]. In this work chains were visualized as deformable droplets. Initially, when a chain arrives 




Figure 11. Results of Barford et al. [186] for irreversible adsorption (black lines) compared to equilibrium 
(gray lines). Scales are arbitrary, (a) Total surface coverage r°° as a function of ft, a dimensionless measure 
of monomer sticking energy, (b) Density profile c as a function of distance from surface, ft = 0.1 and n = 5 
shown on left and right, respectively. Reproduced from ref. [186]. 




Figure 12. (a) Typical configuration of a chain adsorbed during the late stages, according to theory of 
refs. [133,134]. Such chains can adsorb onto free empty sites only (shown as discs) which are separated by 
i scp . In so doing they form loops of s monomers, with as 3 / 5 = ( S cp- (b) Sketch of final irreversibly formed 
layer structure, as predicted in refs. [133, 134], The layer consists of two parts (one chain from each part is 
highlighted): (i) An inner region of flattened chains making as wN contacts per chain, where ui is of order 
unity, (ii) A diffuse outer layer built up from chains each making / jV <g; N contacts with the surface. The 
values of / follow a distribution P(f) ~ / -4 / 5 - Each / value corresponds to a characteristic loop size for a 
given chain, s n cont //. 



at a bare surface, each droplet adsorbs onto a certain maximum cross-sectional area. As available surface 
area for adsorption become scarce, in order for late-arriving chains to adsorb it was assumed droplets deform 
by reducing their cross-sectional areas parallel to the surface to fit into the empty space. In so doing, they 
become more extended into the bulk. Using this model a P(f) distribution was generated whose shape is 
similar to the experimental one of figure HDT bl. 

Motivated by the experiment of refs. [16, 17] on PMMA adsorption, O'Shaughnessy and Vavylonis 
[133, 134] attempted to calculate layer structure at the detail of distributions of loop sizes and contact 
fractions, f2(s) and P(f). In this scaling theory, early arriving chains formed a flattened surface monolayer 
while late arrivers were visualized to adsorb irreversibly onto the few available unconnected empty surface sites 
(see figure IT2T a)). Following ref. [133], if each chain-surface adhesion point involves n cont attached monomers, 
then the surface density of free "supersites" (unoccupied surface patches large enough to accommodate n con t 
monomers) is p supor w Ar bound /ra cont where Ar bound = T^ ound - r bound is the density of available surface 
sites and r b "^ und the asymptotic density of surface-bound monomers. As the surface approaches saturation, 

— 1/2 

the separation between supersites l scp « Psupcr becomes large and the minimum loop size s which can come 
down just connects two free supersites, i.e. as 3 / 5 = Z sop whence s — (n cont /a 2 Arbound) 5 ^ 6 - Chains adsorbing 
at this stage end up as trains of n con t monomers separated by loops of order s units (see figure IT^T a)). The 
bound fraction of those chains is / = dAr bound /dAr ks n cont /s for large s, where Ar is the deviation from 
the asymptotic coverage T°°. Integrating this process up to saturation gives 

a 2 Ar bound = n cont (a 2 Ar/6) 6 , P(f) = AT 4 / 5 (/ « 1) (6) 

where A is a constant of order unity. Equation |((>J| describes a tenuously attached outer layer (small / 
values) of late arriving chains, dressing the flattened inner layer (see figure IT^T b)). Two notable features of 
P(f) are: a peak due to the outer layer at small / and a peak from the early stages centered at / = u. 
These results capture some essential features of the experimental data for PMMA as shown in figure ED 
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Figure 13. (a) Adsorbed polymer mass T versus surface-bound part ri, oun< j as predicted by refs. [133,134]. 
For irreversible physisorption, early arriving chains have a fraction w of their monomers bound to the 
surface and during early stages r\, oun d fa ujT. The early stage prediction for chemisorption is different, 
Fbound ~ r 8 / 3 . The predicted form should be compared with the experimental results shown in figure lToT al. 
(b) Frequency histograms for bound mass fraction, /. Experiment (grey) from figurc lTuT b'). Theory (empty), 
from predicted distribution P(f) ~ / — 4 / 5 with / m i n < / < ui, where values for / m i n = 0.09 and uj = 0.47 
were derived from ref. [16]. 



Interestingly, the predicted scaling forms for the loop distribution and density profile, Q(s) re a~ 2 s~ u ' 5 and 
c(z) ~ z~ 4 / 3 , are the same as those predicted for equilibrium layers (see equations 0) and 0). However, a 
fundamental distinction is that the predicted distribution of bound fractions, figure ITST b). is very different 
to the equilibrium one shown in figure 02b). This reflects the fact that, unlike equilibrium, all chains are no 
longer statistically identical; rather, there are infinitely many families, each parametrised by an / value. 

5. Chemisorption 

Chemisorption is characterized by very small values of the local monomer-surface reaction rate Q. Despite 
the importance of chemical adsorption in many applications, few experimental and theoretical works exist. 
Irreversibility effects are almost always relevant since chemical bond energies are virtually always much larger 
than kT . 

5.1. Theory 

The important feature of chemisorption is that due to the smallness of Q, chemical bond formation requires 
times which are typically longer than the relaxation time of the attached chain. Thus the reaction rate 
k(s) of a monomer s units from a graft point is proportional to the equilibrium probability of finding this 
monomer on the surface, given the current constraints due to all chemical bonds formed at earlier times, 
see figure El The kinetics of single chain chemisorption were first studied by Shaffer et al. [46, 189] who 
focused on PMMA chemisorption onto Al. The form of k(s) corresponding to theta solvent solutions was 
substituted into the master equation describing the evolution of the chain's loop distribution. Electronic 
structure calculations [189] indicated enhanced reaction rates for monomers neighboring a graft point. Due 
to this local cooperative effect, the numerical solution of the master equation led to a process of chain zipping 
onto the surface growing outwards from the first attachments. 

Ponomarev et al. [167] used the bond-fluctuation model to simulate single chain physisorption and 
chemisorption by varying monomer sticking energies and adsorption barrier heights. For chemisorption, in 
the limit of small Q they found r a d s ~ j^o.s±q.2 w j^ n a magnitude greater than the bulk chain relaxation 
time. 

Scaling theory and asymptotic solutions of the chemisorbing chain master equation were developed by 
O'Shaughnessy and Vavylonis [133, 134]. Three general modes of adsorption were found, depending on value 
of the exponent describing the decay of the reaction rate, k(s) ~ s~ 6 . For 9 > 2 the total reaction rate, 
J ds k(s), is dominated by small loop sizes and this was found to lead to a zipping adsorption mode, the chain 
zipping down from the first attachment point. For 9 < 1, since the total reaction rate is dominated by large s 
large loops are more likely to form; this was found to lead to a process of uniform chain collapse, very different 
to zipping. The case 1 < 9 < 2, named "accelerated zipping," was found to be intermediate between zipping 
and collapse; the pure zipping process is effectively short-circuited by the occasional adsorption of large loops 
which serve as nucleation points for further zipping. For non-attractive surfaces exact results were obtained 
for the exponent 9 which can be expressed [134] in terms of other known critical exponents characterizing 
polymer networks [190,191]. For dilute solutions the relation 9 = 1 + v was obtained, where v is Flory's 



Non- Equilibrium in Adsorbed Polymer Layers 



19 




Figure 14. Single chain chemisorption. A chain makes its first attachment with the surface and begins its 
adsorption process. The rate of first monomer attachment is proportional to Q, namely the reaction rate 
given the monomer touches the surface. While for physisorption subsequent monomer adsorption is diffusion- 
controlled, for chemisorption the reaction rate k(s) of the s th monomer away from the first attachment is 
proportional to the equilibrium monomer-surface conditional contact probability. For s -C N, where N is 
length of the tail, k Q/s s where the exponent 9 reflects equilibrium polymer statistics near an interface 
subject to the constraints imposed by earlier reactions. 



exponent [81], giving 9 = 8/5 and 6 = 3/2 for good and theta solvents, respectively. Thus the adsorption 
mode is accelarated zipping in both cases, and the adsorption time was found to be r a d s ~ Q~^N 1 . 

O'Shaughnessy and Vavylonis [133, 134] considered also the kinetics of the build up of a many-chain 
chemisorbed layer from dilute solution. Arguing that the timescale for single chain adsorption was much 
smaller than the time for adsorbing chains to interfere with one another, they predicted a monolayer forms 
in the early stages. This is followed by a diffuse outer layer, attaching at the few available empty surface 
sites, just as for physisorption, figure 1121 The final layer structure was thus of the same form, with two 
chain populations as in figure IT2T b) and with the same loop and contact fraction distributions. However, the 
bound versus total mass relation differs for small coverage, Tbound ~ T 8 ^ 3 as in figure ITST a') . The kinetics of 
layer formation are also very different: in physisorption diffusion-control was predicted to apply at all times, 
r ~ i 1 / 2 , while for chemisorption Ar ~ t~ 1/>5 at long times. 

5.2. Experiment 

In contrast to physisorption, very few experiments have examined the fundamentals of polymer chemisorption. 
In principle however, since timescales are intrinsically much longer, the kinetics might be easier to probe. 
Schlenoff et al. [192] studied PS with functional thiol groups along its backbone chemisorbing on gold from 
THF solution. Surface coverage was measured as a function of time using the quartz crystal microbalance and 
the degree of polymer functionalization, q, was varied. Pure PS did not adsorb; for q = 1% and MW=127,000, 
the asymptotic surface coverage was higher and was reached faster than when q = 5%, MW=304,000. Lenk 
et al. [193] studied thiol-functionalized PMMA chemisorption on gold. Similarly to the results of ref. [192] it 
was found that less polymer adsorbed with increasing q. In every case however the final adsorbed amount 
was higher than the amount of unfunctionalized PMMA physisorbed onto gold. The final layer was found 
to contain unreacted thiols. Tsao et al. [194] studied thiol-functionalized PDMS adsorption onto gold from 
dichloromethane. Similarly to refs. [192, 193] less polymer adsorbed with increasing q. However unlike the 
PMMA system [193] the final layer did not contain unreacted thiols except for the case q = 1. These thiol 
experiments are difficult to interpret at present since thiols not only chemisorb via their S group, but are 
also very powerful physisorbers onto gold [195]. In reality, therefore, a complex mixture of simultaneous 
physisorption and chemisorption is involved. It has also been suggested that they may self-associate in 
solution [192]. 

Cosgrove et al. [196] studied the simultaneous physisorption and chemisorption of poly(methylsiloxane) 
on Al by Fourier-transform infrared specroscopy measuring total adsorbed mass as a function of polymer 
concentration in the bulk. Higher final coverages were found with increasing bulk concentration for a 
fixed adsorption time of one week. Physisorption processes rapidly formed a layer on the surface in less 
than an hour while chemisortion processes continued for days. For this system it was found that less 
mass was adsorbed when, prior to the experiment, the surface was initially exposed to a solution of the 
same polymer at lower concentrations. Qualitative agreement with this experiment was found in a lattice 
monte-carlo model developed in refs. [197, 198]. The simultaneous physisorption and chemisorption of 
poly(hydrogen methylsiloxane) on quartz and porous silica was studied by neutron reflectometry and small 
angle neutron scattering in ref. [199]. Adsorption was allowed to proceed for approximately one day at 
different temperatures. It was found that at higher temperatures the final layer had a more extended density 
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profile and the total coverage was higher than at lower temperatures. 

6. Concentration Effects: Irreversible Adsorption from Melts and Semi-Dilute Solutions 

In this section we review irreversible layer formation from polymer melts and solutions at finite concentration. 
This is a common situation arising in many applications involving reinforcement of polymer-solid or polymer- 
polymer interfaces [2,6,12,200] where, for example, surfaces may be functionalized by various types of coupling 
agents [6,201,202]. In general, interfacial adsorption from melts and concentrated solutions may arise from 
cither physical or chemical interactions. Clearly, the behaviour of chains in the immediate vicinity of a surface 
is relevant and such effects have been probed in experiments studying non-equilibrium polymer dynamics in 
thin melt films in contact with surfaces [177,203,204]. Here we we do not review this large and developing 
reasearch field. 

6.1. Physisorption: The Guiselin Brush 

Theory. Guiselin [55] proposed the following experimental procedure to study irreversible polymer 
adsorption from melts. He considered the situation where a melt is exposed to a surface so attractive to the 
polymer chains that they adsorb instantaneously and irreversibly. This can be realized in practice for very 
strong physisorption (e.g. strong hydrogen bonding). Assuming that chain configurations are equilibrated 
before the melt is brought in contact with the surface, then adsorption results in the irreversible freezing in 
of the equilibrium configurations of the chains which touch the surface. Given the ideal chain statistics in a 
polymer melt [81], the resulting loop distribution was found to be 

fi(») w l/s 3/2 , (s < N 1 ' 2 ) . (7) 

Here f2 is proportional to the probability that a random walk originating at the surface makes its first return 
to the surface after s steps. Guiselin proposed measuring this distribution by washing away unattached chains 
with a good solvent (see figure IT5)l . This results in a swollen and stretched polydisperse brush with a density 
profile [55] 

a 3 c{z) a (a/z) 2/5 , (z < h = aN 5/6 ) , (8) 

where h is layer height. An important difference between Guiselin's layer and equilibrium layers in dilute 
solutions (see equation ijTJl) is in the value of the total surface coverage, T = J dzc(z). For equilibrium 
layers, from equation one obtains T a a~ 2 for e > kT. For the Guiselin brush on the other hand, 
most monomers in the melt within the coil radius aiV 1 / 2 from the surface belong to chains which became 
attached. This leads to much higher surface coverages, T « N x / 2 a~ 2 . If the above melt experiment is 
repeated in a semidilute solution of polymer volume fraction <j>, Guiselin predicted that a fraction of order 
unity of monomers within the coil radius [81] R ~ aTV 1 / 2 ^" 1 / 8 from the surface will become attached. This 
was shown by Marques and Joanny [205] to lead to a surface coverage 

a 2 r a (R/a)4> « TVi/207/8 j (9) 

from which the melt case is recovered by setting cf) — 1 (see also ref. [206]). For distances larger than the 
bulk correlation length of the semidilute solution, z > a0 -3 / 4 , the density profile after removal of unattached 
chains was predicted to decay as z -2 / 5 as in equation (JSJ). 

The structure of Guiselin's brush in contact with a polymer solution has been studied theoretically by 
Aubouy and Raphael [207] and experimentally by Auroy and Auvray [208]. Adhesion properties have been 
reviewed by Leger et al. [11]. 

Numerical Simulations. The kinetics of irreversible physisorption from a semi-dilute solute were studied by 
Monte Carlo simulations by Zajac and Chakrabarti [120] assuming instantaneous and irreversible adsorption 
of monomers upon contact with the surface (values of N employed were 50,100, and 200). It was found 
that for long times the surface coverage reached values consistent with an approximate law T ~ jv°- 38 </> - 38 . 
The density profile of the adsorbed chains within the semidilute solution (i. e. without removing unattached 
chains) appeared to decay more slowly than z~ 2 ^ 5 . The first moment of the profile was found to increase as 
0O.2o±o.O3jyo.8i±o.O5 j^- wag f ounc i that the layer consisted of two regions: a loop-dominated inner region and 
a tail-dominated outer region. 

Jia and Lai [210] used the bond-fluctuation Monte Carlo model to study instantaneous and irreversible 
polymer adsortpion from semi-dilute solutions. Contrary to Zajac and Chakrabarti [120] they found that T 
asymptotes its long time value with a stretched exponential time-dependence. Their results for layer height 
of the swollen layer were consistent with Guiselin's prediction in equation JSJ, though the chain lengths which 
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Figure 15. Irreversibly adsorbed layers from polymer melts. Theoretically predicted structures for 
physisorption and chemisorption are rather similar except at small scales, (a) Hierachichal loop structure 
after adsorption from the melt. For slow chemisorption [209] successive loop size scales are frozen in as their 
critical density is reached. Beyond this their formation would force chain stetching and hence is strongly 
suppressed. After time t the loops bigger than s* ~ l/(Qi) 2 get frozen in. For irreversible physisorption [55] 
by comparison all loop sizes freeze immediately. The end result is the same loop distribution but with s* 
essentially one unit, (b) The experiment proposed by Guiselin [55] entails swelling the adsorbed structure 
shown in (a) with solvent and washing away unattached melt chains. The blob structure of the swollen 
layer in good solvent is shown. It consists of two regions separated at z* ~ s * 3 / 5 corresponding to loop size 
s* [55,209], each region having a different density profile. For irrversible physisorption from a melt there is 
no inner region. 



were used were too short to probe the density profile (N < 50). They observed loop and tail regions similarly 
to ref. [120]. 

Experiment. Experimentally, Cohen-Addad [211] and Cohen-Addad and Dujourdy [212] established that 
large surface coverages can be obtained in adsorption from melts. These studies involved hydrogen bonding 
of PDMS melts of various molecular weights onto fumed silica particles. The particles were dispersed in the 
melt for varying periods of time. After washing away unattached chains with solvent, drying, and weighing, 
r was extracted as a function of N. For large enough times, equation © was satisfied with <j> = \. However 
the asymptotic surface coverage was obtained after a very long time (weeks). These slow kinetics were 
interpreted within a model in which water molecules are initially adsorbed onto silica and prevent PDMS 
monomers from adsorbing; the kinetics was thus controlled by the diffusion of water molecules away from the 
surface [212-214]. Other models have also been suggested [215]. 

Auvray et al. [216,217] used small angle neutron scattering to study physisorption of PDMS onto silica 
through strong hydrogen bonding. Samples were prepared by immersing porous silica in a PDMS melt or 
a semidilute solution of PDMS in dichloromethane for a few days. The samples were then washed with 
dichloromethane and measured. By varying both N and the concentration of the semidilute solution, their 
measurements of surface coverages were in good agreement with equation © . Measurements of the density 
profile for adsorption from a melt showed layer heights consistent with h ~ TV 5 / 6 and evidence was found 
that the density profile could indeed decay as a power law. 

Durning et al. [218] used neutron reflectivity to study physisorption of PMMA from melts onto flat 
surfaces of hydroxylated quartz via hydrogen bonding. The melt was spin-cast onto the quartz plates and 
annealed for different time periods at 165°C. Unbound chains were subsequently washed with benzene after 
cooling and the sample was measured. The total irreversibly adsorbed amount T was found to increase 
with annealing time, reaching a plateau after a few days. Since the process of spin-casting does not 
immediately generate equilibrium chain configurations at the interface, the dependence of T on annealing 
time was intepreted to be due to slow equilibration processes occuring at the interface during annealing. 
For the largest annealing times a scaling T ~ 7V°- 47±0 05 was found, consistent with equation JJjJ. Their 
measurements of the density profile of the swollen layer were consistent with a power law decay, with an 
exponent whose value was about twice the predicted value of equation (JHJ , while measured layer heights were 
smaller, h ~ ArO-72±o.03_ 

Marzolin et al. [219] also used neutron reflectivity to study hydrogen-bonding of PDMS and deuterated 
PDMS (dPDMS) onto silicon wafers. Wafers were incubated in a PDMS melt for 12-24 h at 110°C 
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and unattached chains were rinsed by solvent (deuterated toluene for PDMS and octane for dPDMS). 
Measurements of density profiles for dPDMS showed good agreement with Guiselin's prediction, equation 
(JHJ), but PDMS layers appeared to be less swollen. 

The above experiments raise two important issues which deserve further study: (1) Guiselin's thought 
experiment which envisions the instantaneous freezing-in of equilibrium melt configurations in practice may 
require annealing over large timescales [218] for many experimental realizations. If the Guiselin predictions 
are to apply, this timescale, which depends on these unknown chain dynamics, must be large enough for any 
non-equilibrium configurations generated at the interface during melt deposition to have dissapeared. (2) 
After swelling by solvent, one expects desorption-readsorption events and slow chain movements. This will 
eventually transform equation JSJ into the equilibrium profile of equation 

6.2. Chemisorption: The Slow Guiselin Brush 

Since chemical bond formation requires many surface-monomer collisions, chemisorption from melts can be 
thought of as a slow version of Guiselin's case which involved immediate bond formation. This process 
has been studied theoretically by O'Shaughnessy and Vavylonis [209]. Single chain chemisorption kinetics 
was analyzed in a similar way to that for dilute solutions (section 5.1). The chief difference is the reaction 
exponent 9 determining the surface reaction rate k(s) of the s th monomer from a graft point (see figure ITU) . 
The exponent 9 now represents ideal melt statistics due to screening of excluded volume interactions [81] 
leading to a value 9 = 1/2. Values of 9 less than unity lead to collapse kinetics, during which each loop 
splits into two daughter loops of approximately equal size, very different to zipping mode. Since the lifetime 
of a loop of length s is t s w 1/ J Q S ds'k(s') « l/(Qs 1//2 ), by time t the maximum loop size which could have 
survived is s* — l/(Qt) 2 . Thus at time t a single chain consists roughly of ~ N/s* loops of size s*. 

The single chain collapse kinetics can only be realized for very dilute reactive chains in a bulk of identical 
but unreactive chains. When all chains are functionalized, it was shown in ref. [209] that uniniterrupted single 
chain collapse does not proceed to completion. The reason is that if all l/N 1 ' 2 chains per surface site were 
able to collapse independently, each would generate N/s* loops of size s* by time t, giving a net density 
~ N 1 / 2 / s* . This would exceed the critical level, 1/s* 1 ^ 2 , beyond which adsorbed loops would have to be 
strecthed [220]. In refs. [209] it was assumed that loops are generated up to the critical density and no 
further. This leads to a layer which at any time consists of two parts: (i) an outer layer of frozen large loops 
at the critical density, and (ii) an inner layer of loops which have not yet collapsed. If reactions are stopped 
at tfinai and the layer is swollen with solvent, the predicted layer height is h = aN 5 ^ 6 and the swollen layer 
has two regions with different density profiles c(z): c(z) ~ z*~ 2 ^ 5 for z < z* and c(z) ~ z~ 2//5 for z* < z < h 
where z* = (Qifinai) Note that the outer layer has the same power law profile and height as Guiselin's, 
equation JSJ, but the inner part is different. 

7. Discussion 

The experimental studies of strong physisorbing systems from solution in refs. [16, 17] have identified 
an important layer property, namely the distribution of chain bound fractions, P(f). Features of this 
fundamental layer characteristic may be more accessible to experiment than density profiles. In the appendix 
it was shown that for equilibrium layers in the long chain limit this distribution has a sharp peak at a value of 
/ of order unity. Theoretically, then, the cquilbrium layer consists of tightly bound and statistically identical 
chains, with the exception of a few more weakly attached chains (which however determine layer height). This 
general picture is supported by computer simulations [73, 101] which show that in practice for finite length 
chains the peak of P(f) may be rather broad. The experimental discovery of a qualitatively very different and 
broad P(f) for the strongly physisorbing PMMA-oxidized silicon system [16,17] has demonstrated that such 
non-equilibrium layers comprise many classes of chain, each with its own / value corresponding to its time 
of arrival at the surface. In particular, a large subpopulation consists of weakly bound late arriving chains. 
The existence of a broad P(f) is supported by theories addressing irreversible adsorption [17, 133, 134]. 

A promising approach for future experiments would be to probe P(f) and its dynamics for less sticky 
systems which exhibit measurable aging effects, such as PS on oxidised silicon [128, 129, 135]. It would be 
interesting to test what the initial P(f) is in this case and if it does tend to relax towards equilibrium. During 
this relaxation process, the late arriving, weakly bound chains will progressively become more attached to 
the surface. Since weakly bound chains are presumably more easily displaced by bulk chains, it is possible 
that it is the relaxation process of P(f) which underlies the slowing down of exchange rates with increasing 
aging in many of the displacement experiments (see section 3.3). 
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Another promising future probe of processes in non-equilibrium layers lies in measurements of the 
dynamics of the lateral size of individual chains. As the chains which build up the layer flatten down onto the 
surface, theory has suggested that their initial size is of order the bulk size, R ~ TV 3 / 5 . In equilibrium layers 
however, theoretical arguments [100] support a picture in which most chains have a latteral size R ~ A 1 / 2 
with the exception of a small fraction of weakly attached chains extending far into the bulk. This leads one 
to expect a process of chain shrinkage during layer aging by a factor which may typically be of order unity. 
This change could be measurable with smart chain labeling. 

Recently, Zhao and Granick [221] measured chain diffusion coeffients, D, as a function of surface coverage, 
r. They employed PEO chains with MW=10800 adsorbed onto rather weakly attractive surfaces for which 
they estimated e ea O.bkT. They found that D was enhanced when the surface coverage increased towards 
what seemed to be the chain overlap threshold. But for layers whose coverage was far above this limit, they 
observed a large reduction in D. An interesting question motivated by this experiment is understanding how 
the curve D(T) is modified with increasing surface stickiness. 

The differing structures of the equilibrium and non-equilibrium layers should have important 
consequences for layer mechanical properties and interactions between layers. In the picture of refs. [133, 134] 
the scaling of the density profile is the same for fully formed equilibrium and irreversibly formed layers. Now 
the interactions between two polymer-covered surfaces are believed to depend primarily on the density profiles 
[30, 146, 147]. Thus, when two such surfaces are brought into contact using a surface force apparatus [52], for 
example, this suggests force profiles should be qualitatively similar regardless of whether the polymer layers 
formed reversibly or irreversibly. Since the equilibrium layer chains rearrange themselves relatively rapidly, 
differences may manifest themselves at long times or very high forces when significant chain rearrangements 
can occur [48,51,52], leading to chain bridging and other complex effects. On the other hand adhesion 
properties may be very different in the two cases. Suppose a polymer layer is brought into contact with 
an elastomer or a polymer melt which is subsequently cooled [11]. Even though equilibrium and irreversible 
layers have similar loop distributions, the statistics of surface-chain attachments are very different (see sections 
2.2.2 and 4.2). Consider a loop size s large enough to efficiently couple the surface to the bulk phase. For an 
equilibrium layer, this loop belongs to a chain having order N contacts with the surface; for an irreversibly 
formed layer the number of contacts is of order N/s, since the loop distribution for a given chain is roughly 
monodisperse. This suggests the equilibrium layer may produce a significantly stronger interface, though the 
details are complex and depend on the mode of fracture [12]. 

Finally, the process of polymer chemisorption remains largely unexplored experimentally and many 
possibilities for future work exist. While many chemisorbing systems may in practice involve some degree 
of simultaneous physisorption, there are also certain simplifications since the chemisorption process is very 
slow and often very close to the clean limit of complete irreversibility Future research on the structure 
of chemisorbed layers promises to yield important insights into the many chemisorption-based applications 
involving reinforcement of polymer interfaces and adhesion. 

Acknowledgments 

This work was supported by the National Science Foundation under grant no. DMR-9816374. 

Appendix A. Calculation of Equilibrium Distribution of Bound Fractions, P oq (f) 

In this Appendix we calculate P° q (f) for a chain of length N in an equilibrium layer based on the statistics 
of loops and tails. We will make the approximation that the partition function Z(n\N) of a chain of length 
N making n contacts with the surface (i. e. having 2 tails and n—1 loops) is the product of n + 1 "one step" 
loop and tail weighting factors (see end of appendix for a self-consistency argument). 

Definition of partition functions. Denoting Wt(s\N) and Wi(s\N) the weighting factor for a tail 
and loop of size s, respectively, we have: 

(n-i \ n-i 
N ~t 1 -t 2 -Y J s l )Wt{t 1 \N)W t {t 2 \N)\{w l { Sl \N) . (A.l) 
»=1 / i=l 

The authors of refs. [31, 100] evaluated the scaling forms of the partition functions Zi(s\N) and Z t (s\N) of 
loops and tails in the layer, respectively, which generate the self-similar equilibrium density profile. Thus the 
properties of W* and Wi in equation (|A.1|I must be such that 

oo 

Z t (s\N) = W t {s\N) Z{n\N - s) , 

71=1 



Z(n\N) 



N 



N 



dt- 



n 

i=l 



N 



dsi 
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oo oc pN-s 

Z l (s\N)=W l (s\N)J2Y, ds 1 Z(n\s 1 )Z(m\N - s - s x ) , (A.2) 

n=l rn=l ^° 

where Z, defined similarly to equation (|A.1() . is the partition function for a chain with with just one tail (its 
other end being in contact with the surface). We now make the ansatz that Wi and Wt have the same power 
law scaling as Zi and Zt, respectively, and we will later verify self-consistently that equation (|A.2|I is indeed 
satisfied in the limit N — > oo. Using the results of refs. [31,100] we thus consider Wi(s\N) — /i s Air/i(s\N) 
and Wt(s\N) = fj, s A t r]t(s\N). Here Ai and /x are constants of order unity, A t ~ N^ - "'' 4,11 , and 



m(s\N) 



s- Y - lv (s«s*) 



The functions rji and r\t are normalised to unity, 7 » 1.16 is the susceptibility exponent [81], s is continuous, 
and a small-s cutoff at s of order unity is understood. 

Calculation of Z(n\N). We introduce a new variable M in equation (|A.1J| which allows us to simplify 
it by Laplace transforming M — > E as follows: 

Z(n|A0 = //^A^SivHAO , 

/•M /-A/ »-l I" /-M 1 / ri-l 

E M (n\N) = dh dt 2 TT / ds 
Jo Jo ■■, Jo 



5 [M-tr-ti-^Si ] m (t 1 \N)r k (t 2 \N)l[ m (s i \N) . (A.4) 



i=i 



Since 2m is a convolution over loop and tail lengths, one has 
E E (n\N) = [ Vt (E\N)) 2 [ m (E\N)) n - 1 



e (n-l) lmn(.E|iV)+2hi7fc(E|A0 



PS e (n-l)[7j,(£;|JV)-l]+21iuj t (B|JV) (_E<1), (A.5) 

where we used the fact that for E <C 1, one has 77; (£7 1 AT) — > 1. Because the quantity r]i(E\N) — 1 appearing 
in equation (| A. 5|> is small-scale dependent (due to the fact that the exponent l+2u in equation (|A.3|I is larger 
than 2), it is more natural to evaluate instead the Laplace transform of the distribution 5^(n|iV), where 



Y = M-(n-l)si-2s t , si= dssr)i(s\N) , s t = / ds s rj t (s\N) » s* . (A.6) 

Jo Jo 

Note that the loop size s; is of order unity and depends on the cutoff of rji in equation (| A. 3|> , while s t ~ N 1 ^ 2 ^ . 
As will become clear below (equation JA.lOJl l in the limit M — > oo the probability of negative Y values becomes 
zero. For M 3> 1 one can neglect negative Y and after using equation ljA.51) one has 

E' E (n\N) = 3 E (n\N)e E t n - 1 '> 9t+2Ei ' = e («-i)['7 ! (^|JV)-i+s« i ]+2pnn«(i5|JV)+JJ St ] (£ , <<; ^ > (A ?) 



where -B is now conjugate to V. Laplace transforming equation (|A.3(1 one obtains the following expressions 
for the quantities appearing in equation l|A.7|) : 

Using equation (|A.8|I in equation (|A.7|) one finds that the loop and tail contributions dominate for large and 
small E, respectively: 



\m~ ) e- B{n - 1)E (l«E-i«s*) 
* E (n\N)~{ ^ c(Es , r (£r i >0 (A.9) 



where B and C are positive constants of order unity. In deriving equation (|A.9|) we have considered f — n/N 
of order unity since, self-consistently, the probability of very small / values is very small (see equation ©). 
Laplace inverting one has 

^hao « ^ c , c(«) - { ; f: 5 , (a.i 0) 
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where £(x) is a cut-off function. This distribution has a peak at Y ss n 1 ^ 2l '\ a long tail for large Y, and is 
exponentially small for Y <C re 1 ' ^ 2v '. One thus obtains Z(n|iV) by using Y — ► A" — (n — l)s; — 2s t in equation 
(IA.10|I and substituting into equation (|A.4|I : 

ZfnlAO ~ u N A 2 A n ~ 1 - C ( N -( n -^i- 2 ^ \ , A n) 

Z{n\N)~v A t A t [N _ {n _ m _ 2rSt]1+ 2^[ „i/ (2 „) ) ■ ' ' 

This expression will immediately yield P eq (/) (see below). 

Self-consistency check and the value of A[. We now show that our choice of weighting factors 
in equation l|A.l|) leads to loop and tail partition functions in equation (|A.2|I which are consistent with the 
findings of refs. [31, 100] provided the normalisation constant Ai in the 1-loop distribution function obeys 
certain constraints. 

Consider first Zt in equation l|A.2(l . Repeating the above calculations of equations I|A.1|I - (|A.10J| for the 
case of a chain with one tail, one can show that Z(n\N) = p, N AtA^'~ 1 'EN(n\N), where En has the same 
scaling structure as 2jy. Using the definition of Wt(s\N) one thus has from equation (| A. 2|> : 

oo 

Z t (s\N) = fx N A 2 r, t (s\N)J2^-s(n\N - s)e^ lnAl . (A.12) 

n=l 

Now suppose e NlnA i ^> i ; in which case the exponential favours large values of n. Since Sjy_ s has a sharp 
cut-off for n values of order (N — s)/s~i, thus the sum in equation l|A.12fl is then dominated by n values of 
this order. This would imply Z t (s\N) ~ e -( sinA i)/ s i [ s exponentially suppressed for larger s, inconsistent 
with the tail partition function derived in [31, 100]. The latter has a power law dependence on s (note that 
Zt in equation (|A.2|I involves a sum over configurations of the whole chain and unlike Wt it involves a factor 
p N instead of p s ). Thus, self-consistency requires that the value of A\ is such that e NlnAl is non-diverging 
as N — > oo. 

The case e NlnA ' «C 1 is also inconsistent; then the chain partition function, Z(n\N), would be largest 
for n values much less than N (see equation QA.lljl ). This would imply that the mean bound fraction is much 
less than unity, inconsistent with the equilibrium layer for which a fraction of order unity of the adsorbed 
mass contacts the surface. 

Overall, we conclude e N ln Al is of order unity. In this case the sum in equation i|A.12l) is to leading order 
independent of s and thus Z t (s\N) has the same scaling as i] t (s\N); this is consistent with the results of 
refs. [31,100]. Repeating the above arguments for the case of loops, and considering e NlnA i of order unity, 
one finds Zi(s\N) has the same scaling with respect to s as r/i(s\ N). This is also consistent with refs. [31, 100]. 

The distribution of bound fractions. The probability P(n\N) to have n loops is obtained from 
P{n\N) = Z{n\N)/YZ=i Z(n\N). Since the bound fraction is / = n/N, then using P cq (f)df = P(n\N)dn 
and v w 3/5 one has 

N- 1 / 5 r n 

P c *(f) « Xef NlnA < — — - p (/ max - / )A^/6 ( (A.13) 

(/max J ) ' L J 

where X is a normalisation constant, / max ~ and p is defined in equation (J3J. Since e NlnA i is of order 
unity, from equation i|A.13() one recovers equation of the main text. 

Self-consistency of equation (|A.1|) The assumption of independent loop statistics is consistent with 
the hypothesis of ideal lateral blob statistics at every scale in the adsorbed layer [100]. The approximation 
is also self-consistent in the sense that the probability that a chain makes two overlapping loops of the same 
size s at the same spot on the surface is small. The number of successive loops n needed to make a loop of 
size of order s or larger is n ks uring these steps the chain would have performed a 

random walk on the surface of size n 1 / 2 ps s 3 / 5 . Thus the origin of the next loop of size s is just outside the 
region occupied by the previous loop of the same size. 
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